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Abstract 


A theory for representing a parameter- dependent function as a function 
trajectory is described , along with the theory for determining a piecewise 
analytic fit to the trajectory. An example is given that illustrates the 
application of the theory to generating a smooth surface through a discrete 
set of input cross- section shapes. A simple procedure for smoothing in the 
parameter direction is presented , along with a computed example. Application 
of the theory to aerodynamic surface modeling is demonstrated by applying it 
to a blended unng- fuselage surface. 


Introduction 

Analytic curves, such as cubics, parabolas, and 
least square polynomials, arc often used to fit (pre- 
cisely or approximately) an ordered set of points ob- 
tained by experiment or computation. These fit- 
ting methods can also be extended to generate space 
curves to fit discrete three-space vector locations 
(ref. 1). 

One purpose of the present analysis is to demon- 
strate that with little theoretical difficulty, the meth- 
ods can be extended to vectors in higher dimensional 
spaces, and even to functions. In the latter case, 
the resulting polynomials in function space can be 
used to interpolate solutions of a differential equation 
that depends on a parameter, where these solutions 
have been obtained for a relatively small set of values 
of the parameter. However, the primary application 
that is explored here is that of fitting a smooth sur- 
face through a given set of cross-section shapes. This 
surface is expressed in analytic form, and it is gener- 
ated automatically. Thus, it does not involve time- 
consuming hands-on procedures, like those required 
for fitting surfaces with local patches. 

Symbols 

B function space 

Bj polynomial coefficients, may be 

scalars, vectors, or functions, 
i = 1,2, ... 

F function of two variables 

/, y generic functions 

k number of components of joint 

vector 

M coefficient matrix (eq. (3)) 

rn total number of joint vectors 

P, joint vector, ith joint 


R real number line 

t trajectory parameter 

V matrix whose terms are linear 

combinations of joint vectors or 
functions 

X vector position on trajectory in 

(n + 1) -dimensional space formed 
as cross product of axis t with 
n-dimensional space 

x position vector of point on 

trajectory 

x / matrix obtained from prescribed 

tangent vectors at. joints 

x, y, 2 (^artesian coordinates 

A blending parameter, 0 < A < 1 

£ independent variable 

r normalized trajectory parameter 

(cqs. (5), (6), and (8)) 

Analysis 

Curve Fitting in Three-Dimensional Space 

In general, a parametric space curve is expressed 
in the form 

m 

x(i) = < y{t) > (1) 

1 z(t) I 

where t is the trajectory parameter. The theory for 
fitting such a continuous curve through a given set of 
discrete points in space is discussed in reference 1. 
If the point coordinates are specified at discrete 
values of a parameter, then a continuous function 
of t can be synthesized that includes the specified 
points. However, if only a discrete ordered set of 
points is given without reference to a parameter, 
then it is necessary to parameterize the trajectory 
to define a continuous curve through the points. 
This parameterization is accomplished by basing the 


parameter t on some quantity, such as the cumulative 
chord length between successive points. 

To distinguish the original set of discrete vector 
locations from other points on the continuous trajec- 
tory fit through these points, the original points are 
referred to as joints. (See fig. 1.) This terminology 
is consistent with that of reference 1. 


Blended-Parabola Technique 

This blended-parabola technique can be described 
in the present notation as follows. A parabola Xj is fit 
through three consecutive joints Pj_i, P y -, and P J+1 . 
Its equation is 

xj(t! ) = Bij + B 2 jti + (5) 


Linked Cubics 

If the continuous trajectory is defined by a set of 
cubic curves, each segment has the form 

= B u + B 2 , it + B 3ji t 2 + B 4 yZ 3 (2) 

for the ith interval, where the coefficients B ? in the 
polynomial are vectors. The usual way of determin- 
ing these coefficients is to require the trajectory to 
pass through the joint vectors and match specified 
slope vectors at the joints. If these slope vectors are 
determined by choosing the slopes at the initial and 
final points and requiring second-derivative continu- 
ity at the internal joints, the curve is called a cubic 
spline. Then the following matrix equation for the 
rn tangent vectors (x r -, i = l, m) is obtained (ref. 1, 
ch. 5) 

Mx[, = V,,, (3) 


where 

t - tj_ x 

T\ — J - 

tj + 1 — — 1 


(6) 


Similarly, a parabola x J + J is fit through the points 
Pj, Pj+ 1 , and P j_j_2 with parameter 


T2 = 


t - tj 

tj + 2 ~ tj 


(7) 


Between points P^ and P y+1 the parabolas are 
blended by the formula 


X jJ+l( r ) = (1 - T ) x j[ T l(T)] + tx j + 1 [t 2 (t)] (8) 

where the double subscript on the left-hand side of 
equation (8) indicates blending of the parabolas on 
the jth and j + lst intervals, and 


where M is an m x m tridiagonal matrix of scalars, 
and the column matrix j consists of linear com- 
binations of the components k of the joint vectors 
Pj. Although equation (3) actually represents three 
equations (one for each vector component), it does 
not require three solutions. Inverting the matrix M 
provides the solution for all components 

4,1 = M-*V* (4) 


r = 


t-U 


tj + 1 tj 


(9) 


Equation (8) is actually a cubic in r. No blending 
is performed between the first two points or between 
the last two. 


Trajectories Specified in Terms of a 
Parameter 


In fact, if we specify m points in an n-dimensional 
space, the problem does not increase in complexity, 
as the system is solved by inverting the same m x m 
matrix. 

Determining the tangent vectors at the joints by 
requiring second-derivative continuity may result in 
trajectories with rapid, frequent curvature changes 
(wiggles) and large excursions between the joints. 
This effect can be reduced, with some increase in 
complexity, by the use of rational tension splines 
(ref. 2). If it is permissible to forego the second- 
derivative continuity condition, then the problem can 
be avoided by specifying the tangent vectors at the 
joints by a more conservative method or by using an 
alternative to cubic splines, such as the method of 
blended parabolas (ref. 1, ch. 5). 


Suppose a set of points in two-dimensional space 
is given in terms of a parameter t , which can be 
regarded as time as follows: 

x = f(t) (10a) 

y = g(t) (iob) 

Then, we can construct a Z- axis normal to the 
x-y plane and display the time history of the two- 
space trajectory as a curve in the three-dimensional 
space, as shown by the following equation: 


X(z,y,Z) 


f x(t) \ 

< y(t) 


t ) 


(ii) 


2 


As an example, the equation 



is to be replaced by a function F(£,£). (In a math- 
ematical sense, this step is equivalent to increasing 
the dimension to uncountable infinity.) Thus, in the 
(12a) cubic equation (eq. (2)), the coefficient vectors B* are 
replaced by functions B(£) as follows: 


represents a circle in the x-y plane. However, if we 
take t to be time and display the curve 


< 


( cost ^ 
sin t 
t 


(12b) 


in the three-dimensional (x, ?/, t) space, the helix now 
represents the time history of the circular trajectory 
in the physical two-dimensional (x, y) space. 


Generalization of Theory 


This type of distinction is crucial in the gener- 
alization of the theory, where the parameter is al- 
ways considered as an additional independent vari- 
able. That is, points in an n - dimensional space are 
assumed to be specified in terms of a parameter £, 
but the trajectory is taken in the (n + l)-diinensional 
space with the parameter taken as the (n+ 1) coordi- 
nate variable. The trajectories can be distinguished 
by adopting the terminology shown in the following 
equations: 

f l 


x(0 


X2 (*) I 


(13a) 


l x n (t)) 


' x \ (<) ' 

x 2 (t) 


X(<) = < : } 


x n (t) 

{ t J 



(13b) 


It would even be possible to extend this theory 
to vectors representing infinite sequences (n = oo) 
provided that the distance element series defined by 


d s 2 = + dt2 ( 14 ) 


F(Z, t) = B l (0 + B 2 (Ot + Bz(Ot 2 + B 4 (Ot 3 (15) 

The specific points to be fit by the curve are in the 
function space B , but the trajectory exists in the 
cross-product space B x R, where R is the real line 
specifying the parameter t. This effect is analogous 
to the discussion of equations (10) to (13). 

This kind of cross- product space is utilized in 
path-following techniques (refs. 3 and 4) for the solu- 
tion of a nonlinear differential equation that depends 
on a parameter. In fact, one potential application of 
this theory for curve fitting in function space is to 
interpolate solutions of the equation for values of the 
parameter between those for which the exact solution 
is obtained. This application is not pursued in this 
paper. 

In equation (15), the coefficients are determined 
by matching position and tangent vectors at the ends 


the intervals, by 

the formulas (ref. 1, sec. 

5-3): 

Bi(f) 

= P,(0 



(16a) 

B 2 (0 

= P'(0 



(16b) 

b 4 (0 

= P'i+i(0- 

-2P, + 1 (0 + 2Bi 

i(0 + b 2 (0 

(16c) 

b 3 (0 

= p,+i(0- 

- Bi(0 - b 2 (0 

- b 4 (0 

(16d) 


where P(£) denotes an input (joint) shape, and P'(£) 
denotes the specified slope at the joint with respect 
to r (not £) where 


t - tj 

tj 


(17) 


A further point regarding the computation of the 
coefficient functions should be noted. If cubic splines 
are used, so that the tangent vectors are determined 
by specifying second- derivative continuity with re- 
spect to t at the joints, then a matrix equation has 
to be inverted as in the solution of equation (3). As 
previously discussed for equation (3), no increase in 
complexity occurs as a result of replacing the vectors 
with functions. 


is always convergent. 

However, the purpose of the present analysis is 
not to pursue the development of curve fitting in 
n-dimensional space. Rather, the n-space vector x(£) 


Application to Analytic Surface Modeling 

The problem to be considered is that of fitting a 
smooth surface through a given discrete set of cross- 
section shapes. The shapes are given as functions 
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z i(y) at specified locations x* (i = 1,2 , ...,ra). The 
equations arc simplified considerably by normalizing 
the independent variable ;t/ so that — 1 < y < 1. 


For purposes of illustration, an example was for- 
mulated with five input cross-section shapes having 
analytic descriptions as follows: 

7 ny 

Z\{y) = 0.5 cos -y 

(18a) 

z>(y) = 0.3|1 - i/I -0.4 

(18b) 

~:i (y) = 0.5| sin ny\ 

(18c) 

~i (y) = -0.3|y| -0.25 

(180) 

2 r>(y) = 0.5^/ 1 - y 2 

(18e) 


The cross-section shapes described by these for- 
mulas were' disposed at equal intervals beginning at 
x = 0 and ending at x = 6 and computed at equally 
spaced values of y. These shapes are shown in fig- 
ure 2. The three internal shapes (eqs. (18b) to (18d)) 
have discontinuous derivatives at y — 0. 

Figure 3(a) shows a perspective view of a surface 
obtained by linear interpolation between the input 
cross-section shapes. This surface clearly has slope 
discontinuities in the x-direction. 

Figure 3(b) shows a cubic function fit through the 
input cross sections, with four cross sections shown 
interpolated between each pair of input shapes for 
this particular case. For this surface, the slopes at 
the internal joints were not specified by requiring 
second-derivative continuity at the joints because 
that requirement led to a highly warped surface with 
large excursions between the joints. Consequently, a 
more conservative approach was taken by setting the 
slopes at internal joints equal to the average of the 
left- and right-hand derivates as follows: 


the resulting surface would be bumpy. In this case, 
the input shapes must be preconditioned by replacing 
the data with a smoothed approximation, which can 
be represented analytically in various ways, as for 
example, by smoothing splines (ref. 5). Of course, 
smoothing can be performed only if some deviation 
from the original data is permitted. Once the cross- 
section shapes arc described analytically, whether 
by splines or some other formula, the coefficients 
Bi(€) in equation (15) are analytic functions, and 
consequently, the formula F(£J.) for the surface is 
analytic. 

However, even with such a surface, the variation 
with respect to t may be too extreme for some 
particular application. That is, while the general 
shape of the surface is appropriate, the requirement 
that the surface passes precisely through each of the 
input cross sections may cause too much warping of 
the surface, and a smoother version of the surface is 
required. Such a smoothed version can be generated 
in various ways. One relatively simple procedure is 
described next. 

Assume that the surface is required to coincide 
precisely with the two end sections but that the 
condition of matching the intermediate joint sections 
may be relaxed, If the intermediate shapes are 
completely ignored, the surface that satisfies the end 
conditions with minimum variation in the x-direction 
is a linear blend of the two end shapes: 

z t (y,x) = ( 1 - — ) z\(y) + —zr,(y) (20) 

V T5 / £r> 

At a location x* of one of the joint sections, this 
function has the shape 

z ?i(y) = ( i _ 77^) z \{y) t ~ z T)(v) (21) 

V - F r > / X5 


r/7 

(lx 


^ f z i 4- 1 z i 

2 V •■'•(+! - 


+ 


z i z i. — 1 

x t - x f - +1 


(19) 


With the same input cross sections, the surface 
shown in figure 3(c) was generated by the method 
of blended parabolas. This method also produces 
a smooth fit without large excursions between the 
joints. With a slight increase in complexity, the 
method of rational tension splines (ref. 2) can also 
be used. 


Smoothing 

If the individual input cross-section shapes are 
defined by numerical data rather than analytic func- 
tions, the data may not be smooth, and consequently 


Now this function can be blended with the original 
input shape Zj{y) at x; according to the following 
formula: 


z bi(y) — (1 ^i) z i(y) + \ z ii{y) (22) 

where 0 < A* < 1. Clearly X i = 0 results in 
no deviation from the input section, while A, = 1 
results in total neglect of the input section shape. 
An appropriate value of A, is selected and then the 
function spline procedure is applied with the joint 
shapes z ( {y) replaced with Zfo(y). 

This kind of smoothing is illustrated in figure 4, 
which shows the surfaces that result when the pa- 
rameter A i is varied progressively from 0 to 1. 
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Aerodynamic Surface Example 

Figure 5 shows the steps in transforming the up- 
per surface of a blended wing-fuselage configura- 
tion into its topological equivalent having normalized 
chord lengths. Figure 5(a) is a plan view of the con- 
figuration. Figure 5(b) shows the same configuration 
with the leading edge unswept by subtracting the 
distance from the nose to the leading edge for each 
section, including the fuselage. Figure 5(c) depicts 
the surface with the chord lengths normalized to a 
constant value. 

The function spline procedure was applied to the 
surface of figure 5(c) to interpolate four sections be- 
tween successive joint sections. Then the transforma- 
tions in figure 5 were applied in reverse order. The 
surface grids of the original configuration and of the 
resulting function spline representation are shown in 
figure 6 for purposes of comparison. Figure 7 shows 
a close-up, highly foreshortened, view in the vicinity 
of the wing- fuselage juncture. Figure 8 shows a sim- 
ilar view including both upper and lower surfaces. A 
smooth representation of the surface is obtained. 

Concluding Remarks 

A theory for representing a parameter- dependent 
function as a function trajectory has been presented, 
along with the theory for determining a piecewise an- 
alytic fit to the trajectory. An example is given that 


illustrates the application of the theory to generat- 
ing a smooth surface through a discrete set of input 
cross-section shapes. A simple procedure for smooth- 
ing in the parameter direction was presented, along 
with a computed example. Application of the theory 
to aerodynamic surface modeling was demonstrated 
by applying it to a blended wing-fuselage surface. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
June 2, 1992 
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